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Abstract. In this paper we derive a probabilistic representation of the deter- 
ministic 3-dimensional Navier-Stokes equations based on stochastic Lagrangian 
paths. The particle trajectories obey SDEs driven by a uniform Wiener pro- 
cess; the inviscid Weber formula for the Euler equations of ideal fluids is used 
to recover the velocity field. This method admits a self-contained proof of 
local existence for the nonlinear stochastic system, and can be extended to 
formulate stochastic representations of related hydrodynamic-type equations, 
including viscous Burgers equations and LANS-alpha models. 



1. Introduction 

Formalisms representing solutions of partial differential equations as the expected 
value of functional of stochastic processes date back to Einstein and Feynman in 
physics, Kolmogorov and Kac in mathematics. The theory for linear parabolic 
equations is well developed pTl J6]. The theory for nonlinear partial differential 
equations is not as well developed and it involves either branching processes [19], 
or implicit, fixed point representations. Our representation is of this latter kind, 
'nonlinear in the sense of McKean': the drift term in the stochastic differential 
equation is computed implicitly as the expected value of an expression involving 
the flow it drives. 

Connections between stochastic evolution and the deterministic Navier-Stokes 
equations have been established in seminal work of Chorin [6]. In two dimen- 
sions, the nonlinear equation obeyed by the vorticity has the form a Fokker-Planck 
(forward Kolmogorov) equation. The Biot-Savart law relates the vorticity to the 
velocity in a linear fashion. These facts are used by Chorin to formulate the ran- 
dom vortex method to represent the vorticity of the Navier-Stokes equation using 
random walks and a particle limit. The random vortex method has been proved to 
converge by Goodman [12] and Long [21] , see also [22] . A stochastic representation 
of the Navier-Stokes equations for two dimensional flows using similar ideas but 
without discretization is given in [4]. 

A heuristic representation of the Navier-Stokes equations in three dimensions 
using ideas of random walks and particles limits was done by Peskin [23] . There are 
many examples of modeling approaches using stochastic representations in physical 
situations, for instance [7i |24]. 
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The three dimensional situation is coinpHcated by the fact that there are no 
obvious Fokker-Planck hke equations describing the solutions. LeJan and Sznitman 
[20] used a backward-in-time branching process in Fourier space to express the 
velocity field of a three dimensional viscous fluid as the average of a stochastic 
process. Their approach did not involve a limiting process, and this led to a new 
existence theorem. This was later generalized [2] to a physical space analogue. 

More recently, Busnello, Flandoli and Romito developed a representation for the 
3-dimensional Navier-Stokes equations using noisy flow paths, similar to the ones 
considered here. They used Girsanov transformations to recover the velocity field 
of the fluid from magnetization variables, and generalized their method to work for 
any Fokker-Planck-type equation. Their formulation however does not admit a self 
contained local existence proof, which is possible using our formulation. 

In this work we describe a stochastic representation of deterministic hydrody- 
namic equations that is based on the combination of the active vector formulation 
of the inviscid equation and the shift by a uniform Wiener process. The procedure 
consists thus of two steps. In the first step a Weber formula is used to express 
the velocity field of the inviscid equation in terms of the particle trajectories of the 
inviscid equation without involving time derivatives. The second step is to replace 
classical Lagrangian trajectories by stochastic flows driven by the velocity field. Av- 
eraging the stochastic trajectories produces the solution. It is essential that time 
derivatives do not appear in the expressions to be averaged. 

Perhaps the simplest illustration of this approach is the stochastic formulation 
of the viscous Burgers equation. The inviscid Burgers equation is 

dtu+{u-V)u = 

with initial data 

0) = uo{x). 

The absence of a pressure term allows the velocity to be transported by the fluid 
flow. Thus if X is the fluid flow map, uo X is constant in time and u = uqo (Xf^^) 
gives a Weber formula to recover the velocity from the instantaneous flow map. 
More explicitly, the system 

X = u 

At - X-^ 

U = Uo(A) 

with initial data 

X{a,Q) = a 

is equivalent to the inviscid Burgers equation before the formation of shocks. By 
the notation X^^ ^ we mean the spatial inverse of the map X at time t. 

It can be shown that if we add uniform noise to particle trajectories, then the 
average of the stochastic system satisfies the viscous Burgers equation. (We remark 
that we do not have to restrict our attention to gradient solutions u = V0.) 

Theorem 1.1. LetW be a n- dimensional Wiener process, k^l and uq G (7^+1. 
Let the pair u,X be a solution of the stochastic system 



(1.1) 



dX ^udt + V2iy dW 
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(1.2) 
(1.3) 



A = X 



1 



u = E[mo o A] 



with initial data 



(1.4) 



X{a,Q) 



For boundary conditions, we demand that either u and X — I are spatially periodic, 
or that u and X — I decay at infinity. Then u satisfies the viscous Burgers equations 



with initial data uq. Here E denotes the expected value with respect to the Wiener 
measure and I is the identity function. 

Remark 1.2. The spatial invertibility of X can be seen as follows: An elementary 
computation shows 



showing Xt is locally invertible for all time. Global invertibility now follows since 
our boundary conditions ensure Xt is properly homotopic to the identity map and 
hence has degree 1. We also remark that since the noise Wt is spatially constant, 
the inverse map At is as (spatially) regular as the flow map Xt . We refer the reader 
to pai m E] for the details. 

Remark 1.3. If a solution to the system ( |l.ip -( lL4] ) exists on the time interval 
[0, T], then our proof will show that u satisfies equation ( 11. 5j ) on this time interval. 
Though global existence for ( jl.5) is known, fixed point methods will only yield a 
local existence result for ( jl.lj )-( nL4l ). 

Conversely, given a (global) solution u of ( jl.5j ) which is either spatially periodic 
or decays at infinity, standard theory shows that ( 11.1) has a global solution. Now 
uniqueness of strong solutions for linear parabolic equations and Corollary ( 14. 3j ) 
shows that ( 11.3) is satisfied for all time. We refer the reader to [15] for the details. 

We provide two proofs of the corresponding result for the Navier-Stokes equa- 
tions (Theorem 12.2) , each of which can be adapted to yield the above theorem. 
We present one proof of Theorem 1.1 in Section 51 and leave the details of the sec- 
ond proof to the interested reader. In the next section we describe our stochastic 
formulation of the Navier-Stokes equations. 

2. The stochastic Lagrangian formulation of the incompressible 



In this section we show how the method described in the introduction can be 
used to produce a stochastic formulation of the Navier-Stokes equations. We begin 
by recalling a method of recovering the velocity field from the instantaneous flow 
map of an inviscid fluid. This is the Eulerian-Lagrangian formulation [9] which 
describes the Euler equations as an active vector system. 

Proposition 2.1. Let fc ^ and uq G (7^+1^" divergence free. Then u satisfies 
the incompressible Euler equations 



(1.5) 



dtu + {u ■ V) u — lyAu = 



det(VXi) = cxp / ■ Ugds 
Jo 




Navier-Stokes equations 



(2.1) 
(2.2) 



dtu+ {u-\/)u + \/p = 
V • u = 
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with initial data uq if and only if the pair of functions u, X satisfies the system 

(2.3) X^u 

(2.4) A = X-^ 

(2.5) u = P[(V*A)(mooA)] 

with initial data 

(2.6) X{a,0)^a. 

Here P is the Leray-Hodge projection [8[ [25] on divergence free vector fields. We 
impose either periodic boundary conditions and demand that u and X — I are spa- 
tially periodic, or demand that u and X — I decay sufficiently rapidly at infinity (I 
is the identity map). 

The proof of one direction of this theorem foUows immediately from Lemma 13.31 
A complete proof can be found in [9] . Now to obtain a solution to the viscous system 
we use exactly the same Weber formula ( 12.5) for u, but consider noisy trajectories 
instead of deterministic ones and average the noise out. 

Theorem 2.2. Let v > 0, W be an n-dimensional Wiener process, k ^ 1 and 
Uq e (jk+i,a ^ given deterministic divergence free vector field. Let the pair u, X 
satisfy the stochastic system 

(2.7) dX^udt + V2^dW 

(2.8) A^X^^ 

(2.9) u = EP [(V*A) {uo o A)] 

with initial data 



(2.10) X{a,0)^a. 

We impose boundary conditions by requiring u and X — / are either spatially peri- 
odic, or decay sufficiently at infinity. Then u satisfies the incompressible Navier- 
Stokes equations 

dtu + {u-W)u- v/\u + Vp = 
V • u = 

with initial data Uq. 

Remark 2.3. Remarks 11.21 and 11.31 are also applicable here. See also |T5] . 

Remark 2.4. In the presence of a deterministic external force /, we only need to 
replace uo in equation ( ]2.9) with Lp defined by 

Vt^uo+ I {V'X)f{X„s)ds. 

JQ 

We provide a proof of this in Section \4\, however we remark that the proof of 
Theorem 12.2 given in Section \3\ can also be adapted to yield this. Clearly, if the 
forcing is random and independent of the Wiener process W, then our procedure 
provides a representation of the stochastically forced Navier-Stokes equations. 
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Remark 2.5. The construction above can be modified to provide a stochastic repre- 
sentation of the LANS-alpha (or Camassa-Holm) equations. The inviscid Camassa- 
Holm [5] [rs] equations are 




Lemma 13.3 gives a formula to recover v from the inverse of the flow map. Thus 
we obtain a stochastic representation of the viscous Camassa-Holm equations by 
replacing in p77l)- pTT0l ) with 



The velocity v will now satisfy the viscous equation 

dtv + {u-V)v+ (V*u)u - vAv + Vp = 0. 

We draw attention to the fact that the diffusive term is vAv and not vAu. However, 
we do not derive the relation u — [l~a'^A)~^v; and any other translation-invariant 
filter u — Tv would work as well. 

We provide two independent proofs of Theorem 12. 21 We postpone these proofs 
to sections [3] and \4\ respectively, and devote the remainder of this section to conse- 
quences of this theorem. 

The first consequence we mention is that we have a self contained proof for the 
local existence of the stochastic system ( |2.7) " pTT()|) . 

Theorem 2.6 (Local existence). Let k ^ 1 and uq G (jk+i,a divergence free. 
There exists a time T — T{k,a,L, ||wo||fc+i.a), but independent of viscosity, and a 
pair of functions A, u G C([0, T], C''^^'") such that u and X = I + \ satisfy the 
system pTTll- pTOl l. Further 3U = U{k,a,L, ||uo|U-+i,a) such that for all t € [0,T] 
we have \\u{t)\\k+i.a ^ U. 

Here L > is a given length scale, and C'^'" is the space of {k, a) Holder 
continuous functions which are periodic with period L. The theorem is still true 
if we consider the domain R'^, and impose decay at infinity boundary conditions 
instead. We do not present the proof in this paper, but refer the reader to [l4i 15] . 

We remark again that our estimates, and existence time are independent of 
viscosity. The theorem and proof also work when the viscosity v = 0, thus we have 
a proof that gives us local existence of both the Euler and Navier-Stokes equations. 

The nature of our formulation causes most identities for the Euler equations (in 
the Eulerian-Lagrangian form) to be valid in the above stochastic formulation after 
averaging. We begin by presenting identities for the vorticity. 

Proposition 2.7. Let lo — V x u be the vorticity, and ujq — V x uq be the initial 
vorticity. Then 




dtV +{u-W)v+ (V*u)t; + Wp^O 
u = (1 - a'^Ay^v 
V • w = 



(2.11) 
(2.12) 



u = EP [(V*A) uo o A] 
u = (1 - a^Ay^v. 



(2.13) w = E[((VX)wo)o A]. 

// the flow is two dimensional then the above formula reduces to 

(2.14) tj = E[woo^]- 
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Remark 2.8. In the presence of an external force /, we have to replace coo in equa- 
tions pJ3l l and pJlj ) with w defined by 

zut^ujo+ f {VX,y^giXs,s)ds 
Jo 

where g — \7 x f . For two dimensional flows this reduces to 

= cjo + / g{Xs,s)ds. 
Jo 

We draw attention to the fact that these are exactly the same as the expressions in 
the inviscid case. 

We can prove this proposition in two ways: The first method is to directly 
differentiate ( \2.9\ ). and use the fact that P vanishes on gradients. We leave the 
details of this to the interested reader. The other method (presented in Section 3| 
is to use Proposition 14.21 and the generalized Ito formula to show that ui satisfies 
the vorticity equation 

(2.15) + (u • V) w - i^Aw = (w • V)u. 

The two dimensional vorticity equation does not have the stretching term {lo ■ 
V)u. With our formulation the two dimensional vorticity equation is immediately 
obtained by observing that the third component of the flow map satisfies X^la, t) = 
as- 

Although it is evident from equations ( j2.9j l and ( ]2.13j ), we explicitly point out 
that the source of growth in the velocity and vorticity fields arises from the gradient 
of the noisy flow map X. The Beale-Kato-Majda [1] criterion guarantees if the 
vorticity to stays bounded, then no blow up can occur in the Euler equations. In 
the case of the Navier-Stokes equations, well known criteria for regularity exist 
and they can be translated in criteria for the average of the stochastic flow map. 
As is well known, the vorticity of a two dimensional fluid stays bounded, which 
immediately follows from equation ( |2.14j ). 

Finally we mention the conservation of circulation. For this we need to consider 
a stochastic velocity u deflned by 

(2.16) u = P [(V*A) (uo o A)] . 

We remark that u is defined for all realizations of the Wiener process, because the 
diffusion matrix in equation ^.7) is spatially constant. This can be seen from the 
techniques used by LeBris, Lions [18] or those used in Section[3l Notice immediately 
that u = Eu and uq — uq. The circulation of the stochastic velocity is conserved 
by the stochastic flow. 

Proposition 2.9. IfT is a closed curve in space, then 

u ■ dr — d) Uq ■ dr. 

x(r) Jr 

Proof. By deflnition of P, there exists a function q so that 

(V*A)(uooA)+Vg 

=^ {V^X){uoX) = ua + ViqoX). 
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Hence 



/ u-dr^f {uo X oT) ■ {WX\^r')dt 
Jx(r) Jo 

= [ (v'x|j,)(uoXor)-r'di 



(uo + V((7 o X)) ■ dr = f uq ■ dr. 

We remark that the above proof is exactly the same as a proof showing circulation 
is conserved in inviscid flows. □ 

The rest of this paper is devoted to proving the results stated in Sections Q] and 
[21 We provide two independent proofs of Theorem 2.21 however only adapt the 
second proof to yield Theorems l.lj 2.21 and Proposition 2.71 

3. Proof of the stochastic representation using pointwise solutions 

In this section we prove Theorem 12.2 by constructing pointwise (in the proba- 
bility space) solutions to the SDE ( j2.7j) . This idea has been used by LeBris and 
Lions in [18] using a generalization of the W^'^ theory. In our context however, the 
velocity u is spatially regular enough for us to explicitly construct the pointwise 
solution without appealing to the generalized W^'^ theory. 

We remark that the proof given here will also prove Theorem 1.11 and leave the 
details to the interested reader. We begin with a few preliminaries. 

Definition 3.1. Given a (divergence free) velocity u, we define the operator V by 

Vv = df V + (u • V) u 

Lemma 3.2. The commutator [D, V] is given by 

[V, V]/ = P(V/) - V(I?/) = -V*uV/ 

Proof. By definition, 

[I?,V]/-P(V/)-V(P/) 

= (u-V)V/~V[(u.V)/] 

= (u- V)V/- (u- V)V/- V*uV/ □ 
Lemma 3.3. Given a divergence-free velocity u, let X and A he defined by 

X = u{X) 
X{a,0) = a 
A^X-^ 

We define v by the evolution equation 

Vv = F. 

with initial data vq. If w is defined by 

w = P [(V*A) v] , 
then the evolution of w is given by the system 

dtw + {u ■ V)w + (V'u)u; + Vp = (V*A)r 
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V • W = 

Wo = Pwo 

Proof. By definition of the Leray-Hodge projection, there exists a function p such 
that 

w = V -Vp 
= v.i\IAi - Vp 
=^ Pu) = {Vvi)VAi + UiPVA^ - PVp 

= r,VA, - t;,(V'u)VA, - Wp + (V'-u)Vp 

= (vM)r - (v*u)(wiVA, + Vp) - vvp 

= (v'A)r - (v*u)w - vvp. □ 

Corollary 3.4. // u, X, A are as above, and we define w by 

w^P [(VU) uqoA]. 

Then w evolves according to 

Vw + {V*u)w + Vp = 

V • w = 
w{x, 0) = ^0(2^) 

Proof. The proof follows from Lemma 13.31 by setting v = uq o A and F = 0. □ 

We now return to the proof of Theorem 12. 21 

Proof of Theorem \2.2[ For simplicity and without loss of generality we take — \. 
Let (ri, P) be a probability space and W : [0, 00) x ^ a three dimensional 
Wiener process. Define u'^ and by 

u'^ix.t) = u{x + Wt{ijj),t) 

and 

^^(0,0) = a 

Although is not Lipschitz in time, it is certainly uniformly (in time) Lipschitz 
in space. Thus the regular Picard iteration will produce solutions of this equations. 
Finally notice that the map X defined by 

solves the SDE ( 12771 ). 

Let B'^ be the spatial inverse of . Notice that 

X{B'^{x-Wt,t),t) = Y'^{B'^{x-Wt,t),t) + Wt 
= X 

and hence 

(3.1) A^Tw.B'^ 
where Tx is the translation operator defined by 

(3.2) rj{y)^f{y-x) 
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We define u>" by 

w"" =P [(V*B") uo o B"] . 
By Lemma |3.3i the function w'^ evolves according to 

(3.3) dtw'^ + {u'^ ■ W)w'^ + (V*w'^)w + Wq'^ = 

(3.4) V • = 

(3.5) w'^ix^O) ^uo{x). 
Now using equation ( ]2.9j ) we have 

u = EP [(V*A) Uo o A] 
= EP [{V'tw,B^)uootw,B'^] 
= EP [tw, ((V*i?") uo o i?-)] 
= ErvF.P [(V^B"^) uooB'^] 

(3.6) = Etw.w'^ 

Our assumption mq G (^fe+i^a along with Theorem 12.61 guarantee that w'^ is 
spatially regular enough to apply the generalized Ito formula \17] to w"(x —Wt,t), 
and we have 



w'^ix -Wt,t) -uo{x) = / w'^{x-Ws,ds) 



I X — Ws , s 



Ai 



I a; — Ws ; s 



ds+( djw'^ix -Ws,ds) , Xj -Wj! 



Notice that the process is in time (since the time derivative is given by 
equation (3.3) ), and hence bounded variation. Thus the joint quadratic variation 
term vanishes. Taking expected values and using ( j3.6[ l we conclude 



(3.7) u{x,t) - ua{x) w'^ {x - Ws, ds) + ^ Au{x,s)ds 

Jo Jo 

Using equation ( 13. 3p and the definition of the Ito integral we have 

e/ w^{x-Ws,ds)=E [ dtw'-l_^ ds 
Jo Jo 

= E f [{u- ■ V)w- + (V'u")u;- + Vq-] ^ ds 
Jo 

E I \{u(x,s)-V)w^l^^^^^ + {V'u{x,s))w'^l^^^^^- 
- Vq' 



I X — Ws ■ s 



ds 



(u(x, s) ■ V) u\^ ^ + {W'u{x, t)) ^ + 



VEq'^{x-Ws,s] 



ds 



(3.8) 

where q' is defined by 



[u{x,s)-V)u\+^q'\ 



q' = ^W\u\^ +ETw,q^ 



ds 
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Using equations ( ]3.8j ) in ( ]3.7j ) (along with the observation that the joint quadratic 
variation term is 0) we obtain 

u{x, t) — Ui){x) ~ — I [(u • V)m + ^ Au + Vg] ^ ^ ds 
Jo 

Equations ( 13. 4j ) and ( 13. 6p show that u is divergence free, concluding the proof. □ 

4. Proof of the stochastic representation using the Ito formula 

In this section we provide a proof of Theorems l.lj 12.2 and Proposition 12.71 
directly using the generalized Ito formula. At the end of the section we provide a 
brief contrast with the diffusive Lagrangian formulation [10] . 

Before beginning our computations we remark that our assumptions fc 1, 
uq £ (jk+i,a g^j^j Theorem 12.61 guarantee and that the processes A and X are 
spatially regular enough to apply the generalized Ito formula pT]. Further, as shown 
in p3] the displacements X — I and A — I are spatially C'^+i'" (thus bounded in 
periodic domains) and hence the Ito integrals (and expectations) that arise here are 
all well defined. This said, we assume subsequently that all processes are spatially 
regular enough for our computations to be valid. 

We begin with a lemma leading up to computing the Ito derivative of A. 

Lemma 4.1. There exists a process B of bounded variation such that almost surely 
(4.1) At^Bt~V2^ f VA, dWs 

JO 

Proof. We apply the generalized Ito formula to ^ o X to obtain (almost surely) 
= A{X,,ds) + VA|^^_^ dX, + i df^Mx^J (^^ ' ^. 
^ d,A{Xs,ds), xi~xi, 

ds 



^31 s 



(4.2) = / A{Xs,ds)+ f 
Jt' Jt' 



Vv4L, u + vl\A\^ 

I As ,s I As .s 



2v I VA\^ ^dWs 



^ d,A{X,,ds), Xl-Xl, 



Notice that the second and fourth terms on the right are of bounded variation. 
Since the above equality holds for all t',t, and Xg is a homeomorphism, the lemma 
follows. □ 

We remark that the term in ( 14. 2j ) has a positive sign, which is the anti- 
diffusive sign forward in time. The quadratic variation term in (4.2j ) brings in twice 
the negative Laplacian, thus correcting this sign, and we end up with a dissipative 
SPDE for A as we should. 

Proposition 4.2. The process A satisfies the stochastic partial differential equation 

(4.3) dAt + {u-V)Adt- lyAA dt + y^VA dWt = 

Proof. We apply the generalized Ito formula to AoX as in Lemma |4.1i Notice that 
differentiating equation (j4.1j) , we find the martingale part of diA. Since the joint 
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quadratic variation term in equation ( ]4.2j ) depends only on the martingale part of 
diA, we can now compute it explicitly. 

\7[d,A]^ ^dWs,Wl-Wl, 



J^^ d,A{X,, ds) , XI - Xl^ = -2v l^j^^ 



2v I AA\ 



ds 



Substituting this in equation ( 14. 2) , the proposition follows. 



□ 



Corollary 4.3. Let d be spatially , and differentiable in time. Then the process 
V — i!) o A satisfies the stochastic PDE 



(4.4) 



dvt + {u ■ V)vt dt ~ lyAvt dt + \/2i^Vut dWt = dfi^lj^^ dt 



Remark 4.4. Consider the case when -d is independent of time. Now when v — 0, 
we know that equation ( |4.3| ) is equivalent to the fact that A is constant along 
particle trajectories. Thus any function of A, in particular u, is also constant 
along trajectores and hence automatically satisfies equation HA} . When ^ 0, 
Proposition 14.21 and Corolary 14.31 make the same assertion for noisy trajectories. 

Proof. The corollary follows directly from Proposition 14.21 and the generalized Ito 
formula: 



dvt = 



dt + Vi?| dAt + ^d^jd\^d{A, , Aj) 



^3 It 



yd\^^{yAt)ut + vy'd\^AAt + i^df^^l^dkAldkAi 
{ut ■ V)vt + vAvt dt - 



dt- 



2vVvt dWt 



□ 



The results from sections Q] and now follow as direct consequences of Proposition 
ajand Corollary 331 

Proof of Theorem ] The proof of this theorem follows by setting i!} = uq, inte- 
grating ( 14. 4j ) and taking expected values. □ 

Proof of Theorem \2.2l We provide the proof in the presence of an external force 
/, as stated in Remark 12. 41 Let v — (p o A, w — (V'yl)i; and we compute the Ito 
derivative of w using Ito's formula: 

dwi — {diA) ■ dv + d{diA) ■ v + d {diAj , Vj) 

= diA ■ [{-u ■ V)v + vAv + (V*X) 1^/] dt - v^d.A ■ (Vw dW) + 

+ v-[- {(d^u) • V) A - (u • V)d^A + vAd^A] dt - V2^v ■ (yd.AdW) 

+ 2vdl^AjdkVj dt. 

Making use of the identities 

[u ■ \I)wi = d^A ■ [{u ■ V)w] + [{u ■ V)5,A] • v 

Awi = diA ■ Av + Ad,A ■ v + 2dkiAjdkVj 

diUkWk = V ■ [{diu ■ V) A] 

d,A-[{V'X)\j]=f 
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and 

dkWi = diAjdkVj + VjdkiAj 

we conclude 

(4.5) dw = [-(w • V)u; + i^Aw - {V^u)w + f] dt - V2i^Vw dW. 

Now from equation ( 12. 9j ) we have (almost surely) 
u = Ew + Vq 

=^ It - Mo = E f [-(li • V)w + vAw - (V*?i)u) + f]+Vq 
Jq 

= I [-{u • V)(m - Vq) + iyA{u - Vq) - (V'u)(u - Vg) + /] + Vg 
Jo 

= /" [-(u- V)u + i/Au + /]+Vp 
Jo 



where 



P^q- f [{u-V)q-l^Aq + ^\uf] 
Jo 

Differentiating immediately yields the theorem. □ 



Proof of Proposition \2.7\ Notice first that WX is differentiable in time. We set 
■d — {'S/X)u!o, u} — d o A and apply Corollary |4. 31 to obtain 

dii! + {u- V)w dt - vAlo dt + VtI/VCj dW = VdtX\^UQ{A) dt 

= (V«)(VX)|^L^o(^)di 
= (Vm)cj dt 

Integrating and taking expected values shows that = Eo) satisfies the vorticity 
equation fl2.15j ) with initial data wo- The proposition follows now follows from the 
uniqueness of strong solutions. □ 

4.1. A comparison vifith the diffusive Lagrangian formulation. The com- 
putations above illustrate the connection between the stochastic Lagrangian for- 
mulation ( j2.7j )- p^.l0| ) presented here, and the deterministic diffusive Lagrangian 
formulation. We briefly discuss this below. In [1(J| the Navier-Stokes equations 
were shown to be equivalent to the system 

dtA + (m • V) A - vAA = 

u = V\[V^A)v\ 

dtvp -I- (u • V) w/3 — vAvp — 2iyC^j pf^j^i 

CI. = {VA)-,ldud,A, 

with initial data 

l(x,0) = 

v{x, 0) — Uq{x). 

We see first that A = 'EA. The commutator coefficients C"- in the evolution of 
V compensate for the first order terms in A((V'A)ti). In the stochastic formula- 
tion these arises naturally from the generalized Ito formula as the joint quadratic 
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variation term {diAj,Vj). More explicitly, the equations 

and d{diAj , Vj) = 2i'dljAidkVjdt 

illustrate the connection between the two representations. 
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